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1.  INTRODUCTION 

In  this  paper  a  method  for  determining  radar  emitter  pulse  repetition  interval 
(PRI)  in  a  dense  environment  is  presented.  The  PRI  estimates  can  be  used  for 
emitter  identification.  The  method  is  based  on  statistical  techniques. 

This  work  is  done  wider  task  DST  76/004. 

The  aim  of  an  emitter  classification  system  is  to  identify  possible  threats  by 
analysing  intercepted  signals.  The  importance  of  emitter  identification  in  an 
EW  system  can  be  judged  from  the  fact  that  it  is  a  necessary  first  step  in 
implementing  countermeasures  against  threat  emitters. 

To  deal  with  the  complex  signal  environments  likely  to  be  encountered  one 
emitter  parameter  alone  cannot  be  relied  on  to  achieve  positive  identification. 
Thus,  in  this  paper  it  is  assumed  that  some  analysis  of  the  environment  using 
such  parameters  as  radio  frequency  and  direction  of  arrival  will  be  carried  out 
in  conjunction  with  PRI  analysis. 

The  electromagnetic  environment  can  be  divided  into  two  classes  of  emitters, 
namely,  those  emitters  whose  intercepted  pulse  amplitude  is  heavily  modulated 
(such  as  search  radars)  and  those  emitters  whose  intercepted  pulse  amplitude  has 
a  modest  amount  of  or  no  modulation  (such  as  tracking  radars) .  The  range  of  PRI 
to  be  encountered  will  tend  in  practice  to  be  from  about  102  Ms  to  104  Ms. 

Search  radars  tend  to  have  longer  PRI  while  tracking  radars  usually  have 
relatively  short  PRI. 

Standard  methods  for  emitter  recognition  using  pulse  repetition  interval 
include  the  pulse  train  sorter  and  the  time  of  arrival  method.  Experience  has 
revealed  certain  limitations  for  these  methods.  The  main  limitations  are  in 
handling  a  variable  PRI  and  in  extending  their  use  to  a  dense  environment. 

It  is  important  to  be  able  to  recognize  emitters  with  variable  pulse 
repetition  interval.  Variable  PRI  is  used  in  radar  emitters  to  overcome  the 
blind  speed  problem  with  moving  target  indication  (MTI)  and  as  an  antijamming 
technique.  In  the  MTI  application  the  PRI  may  vary  typically  ±  10%  whilst  in 
the  antijamming  application  ±  50%  has  been  considered. 

The  recognition  methods  in  this  paper  can  cope  with  variable  PRI  and  offer  a 
simple,  systematic  approach  to  analysing  the  dense  environment  problem. 

A  shortened  version  of  this  paper  appeared  in  reference  1. 

The  method  proposed  in  this  paper  will  be  sketched  briefly  now.  In  order  to 
monitor  the  emitter  environment,  data  are  collected  over  a  series  of  consecutive 
observation  intervals.  The  time  between  changes  to  the  environment  of  emitters 
is  assumed  to  be  long  compared  with  an  observation  interval.  Each  observation 
interval  is  divided  up  into  gate  intervals  and,  at  the  end  of  each  observation 
interval,  the  relative  frequency  of  the  various  possible  number  of  pulses  which 
can  occur  during  a  gate  interval  is  recorded.  It  turns  out  that  the  possible 
number  of  pulses  which  can  occur  is  of  the  order  of  the  number  of  emitters  present. 
Using  the  recorded  relative  frequencies  it  is  possible  to  deduce  the  types  of 
PRI  (whether  fixed  or  M-position  stagger,  say)  and  estimate  the  PRI  parameters 
(for  a  2-position  stagger  the  two  PRI  used)  when  more  than  one  emitter  is  present. 
The  data  processing  to  accomplish  this  involves  only  a  few  simple  rational 
operations.  The  method  adopted  is  adaptive  and  uses  a  certain  algorithm  to 
select  the  length  of  the  gate  interval.  Changes  which  occur  to  the  environment 
are  analysed  using  the  relative  frequencies  recorded  at  the  end  of  each 
observation  interval.  In  this  way  the  PRI  parameters  of  the  emitter (s)  which 
started  or  stopped  being  intercepted  are  determined.  This  adaptive  approach 
yields  more  accurate  parameter  estimates  than  is  possible  with  a  non-adaptive 
approach. 

An  example  of  the  achievable  accuracy  is  as  follows.  With  up  to  five  non¬ 
scanning  emitters  present  and  one-second  observation  intervals,  the  r.m.s. 
percentage  error  in  the  estimates  of  fixed  PRIs  is  less  than  0.4%.  The  main 
advantage  of  the  method  is  that  only  a  small  amount  of  data  has  to  be  manipulated. 

The  reader  who  is  not  interested  in  the  mathematics  should  go  directly  to 
Sections  9  and  10  for  a  discussion  of  possible  applications  and  further 
developments . 
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2.  OBSERVATION  PROCEDURE  AND  DEFINITIONS 

In  this  section  a  readily  iraplementable  observation  procedure  is  introduced 
for  the  determination  of  radar  emitter  PRI  from  intercepted  emissions.  Suppose 
that  the  IF  output  of  a  typical  electronic  surveillance  receiver  is  examined 
during  the  time  interval  from  s0  to  s^.  This  interval,  which  will  be  called 

the  observation  interval,  is  partitioned  as  s0 ,  si  ,  . sT  sT  where 

st  -  s  =  t  for  all  t.  The  quantity  r  will  be  called  the  gate  time  and  the 

subintervals  of  length  r  will  be  called  the  gate  intervals.  Note  that 

Tr  =  s^  -  so .  The  observation  procedure  is  to  determine  the  number  of  pulses 

which  occur  during  each  gate  interval  and  set  n^  equal  to  the  number  of  pulses 

which  occur  during  gate  interval  t.  Thus,  a  sequence  of  numbers  ni ,  m ,  ...,  n^, 

is  obtained.  Suppose  that  during  the  observation  interval  the  maximum  and 

minimum  number  of  pulses  detected  in  any  gate  interval  is  N  and  N  . 

*  max  mm 

respectively.  A  set  of  quantities,  x  ,  is  derived  from  the  observations,  n  , 

w  C 

as  defined  by 


x 


t 


n^ 


N 


mm 


(1) 


The  PRI  of  an  emitter  is  defined  as  the  time  interval  between  successive 
pulses.  The  sequence  of  PRI,  yx ,  y2 ,  . . . ,  y^  is  assumed  to  be  a  realization  of 

the  stochastic  process,  {  Y  ,  n  =  1,  2,  . ..,  N],  Four  types  of  emitter  PRI  will 

be  considered  in  this  paper.  Each  type  can  be  characterized  by  its  PRI 
probability  density  function,  pfy^),  as  follows: 

(a)  Fixed  PRI  emitters  with 


p(yn)  =  5(yn  -  z) 


(2) 


so  that  y  =  z  for  all  n. 
n 

(b)  M-position  staggered  PRI  emitters  with 


M 


p(yj  =  M"1  \  «(y  -  Zj 


n 


n  m 


(3) 


m=l 


The  PRI  sequence  is  cyclic.  For  example,  if  M  =  2  then  y2n  =  zi  and 

y2n+l  =  7,2  for  a11  n' 

(c)  Uniform  PRI  emitters  with 


P(yn)  = 


(b  -  a)~  ,  if  a  <  y^  <  b 
0  ,  otherwise 


(4) 


The  PRI  sequence  is  independently  and  identically  distributed. 
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(d)  Sinusoidal  PRI  emitters  with 

P  f  hy  -  a  -  b 

'  Cb  -  aWi  - 

0  ,  otherwise 

The  PRI  sequence  is  independently  and  identically  distributed. 

In  the  above  5(.)  denotes  the  Dirac  delta  function.  When  referring  to  a 

dense  environment  of  emitters,  to  denote  the  k  emitter  a  superscript  k  will  be 
added  to  algebraic  symbols.  For  the  sake  of  simpler  notation  it  will  be  omitted 
where  no  ambiguity  arises. 

Suppose  that  Q  +  1  possible  states  for  can  occur.  The  relative  frequencies 

of  occurrence  of  the  states  of  the  process,  P  ,  q  =  0,  1,  . . . ,  Q,  during  the 

M. 

observation  interval  can  be  calculated  via 


P  =  N  /T 

q  q 


T 


t=l 


N 


t 


q 


1,  if  x  =  q 
0,  if  otherwise 


3.  ASSUMPTIONS 


(6) 

(7) 


(8) 


The  assumptions  made  throughout  this  paper  about  the  radar  emitters  and  the 
observation  procedure  are  stated  in  this  section. 

A  mathematical  model  will  be  developed  for  the  observed  data.  The  proposed 
observation  procedure  generates  pseudo  random  sequences  of  numbers.  The 
equations  derived  for  the  mathematical  model  are  based  on  certain  assumptions 
concerning  these  sequences.  The  validity  of  these  assumptions  should  be 
verified  either  by  theoretical  analysis  or  by  empirical  tests.  (See 
references  2  and  3  for  a  discussion  of  standard  random  number  generators) . 
Theoretical  tests  would  give  much  more  understanding  of  the  method  than  the 
results  of  empirical  tests.  In  this  paper  conditions  concerning  the  PRI  of 
the  emitters  and  the  observation  procedure  will  be  given  which  are  believed  to 
be  sufficient  for  the  mathematical  model  to  adequately  represent  the  observed 
data.  The  results  of  simulations  of  the  proposed  PRI  estimation  technique  which 
have  been  done  justify  this  belief.  Theoretical  proofs  of  the  results  will  be 
given  in  Appendix  I . 

The  basic  assumptions  made  are  as  follows: 

(a)  The  pulse  trains  of  the  emitters  are  unsynchronized. 

(b)  The  pulse  trains  of  the  emitters  are  assumed  to  be  unsynchronized  with  the 
gate  interval  partition  of  the  observation  procedure. 

(c)  If  overlapping  pulses  occur  it  is  possible  to  detect  how  many  pulses 
occurred. 

(d)  The  environment  does  not  change  during  the  observation  interval. 
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(e)  All  pulses  transmitted  are  detected  -  including  those  of  small  amplitude 
in  antenna  side  lobes. 

These  assumptions  are  not  overly  restrictive.  For  emitters  with  fixed  PRI 
assumptions  3(a)  and  3(b)  are  tantamount  to  assuming  that  the  PRIs  of  the 
radar  emitters  are  incommensurable  (i.e.,  not  divisible  without  remainder  by  the 
same  quantity)  and  that  the  PRIs  of  the  radar  emitters  and  r  are  incommensurable. 
In  practice,  radar  emitters  with  commensurable  PRIs  will  be  quite  rare  and  can 
possibly  be  treated  as  a  single  emitter  with  variable  PRI  if  present  together. 
Assumption  3(b)  can  always  be  satisfied  by  suitable  choice  of  r.  It  is 
believed  that  assumptions  3(c),  (d)  and  (e)  are  realistic  in  that  for  practical 
application  the  techniques  can  be  modified  to  deal  with  problems  which  will 
occur.  This  will  be  discussed  further  in  Sections  7  and  9. 


4.  EMITTERS  WITH  FIXED  PRI 


4.1  Single  emitter 

If  a  single  emitter  with  PRI  z  is  present  and  the  gate  time  is  selected 
so  that 


Wz  <t  <  (W  +  1)  z  (9) 

where  W  is  a  nonnegative  integer  then  during  any  gate  interval  either  W  or 
W  +  1  pulses  will  be  observed.  Thus,  the  possible  states  for  x  are  0  and  1 

and  occur  with  probabilities  denoted  by  po  and  p!  =  1  -  p0  respectively. 

An  expression  for  the  probability,  p0 ,  in  terms  of  z,  W  and  r  is  now 
derived.  Let  0^  be  the  time  measured  from  the  start  of  gate  interval  t  to 

where  the  next  pulse  occurs.  Note  that  0  <0^  <  z.  Furthermore,  if 

assumption  3(b)  is  invoked  then  8  may  be  considered  to  be  a  sample  value 

of  a  random  variable  uniformly  distributed  on  the  interval  0  to  z.  The 
proof  will  be  given  in  Appendix  I.  In  figure  1  the  number  of  pulses 
occurring  during  gate  interval  t  is  shown  as  6  varies  from  0  to  z  for  two 

cases  W  =  0  and  3.  In  general, 

x.  =  1,  if  0  <0  <t  -  Wz 

t  ^ 

Xj.  =  0,  if  r  -  Wz  <0^  <  z  (10) 

Since  0  is  uniformly  distributed  on  0  to  z  the  expression  for  po  is 

po  =  [  (W  +  l)z  -  t]  /z  (11) 

Simple  algebraic  manipulation  yields 

z  =  r/{(W  +  1)  -  Po]  (12) 


These  results  are  summarized  in  Table  1. 
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4.2  Dense  environment 

Suppose  that  K  fixed  PRI  emitters  are  present.  The  observation  x  is 
related  to  the  observations  for  the  individual  emitters  comprising  the 
environment,  x^,  via 


K 


k=l 


(13) 


The  possible  values  of  x.  are  the  set  j  0,  1,  ...,  K]  and  K  =  N  -  N  .  . 
r  t  ‘  J  max  mm 

The  probability  that  state  j  occurs  in  a  dense  environment  of  K  emitters  is 

denoted  by  P.  In  order  to  derive  an  expression  for  P  in  terms  of 

J >  K  J 

k  k 

the  quantities,  p0 ,  the  sequences  of  states  x  are  assumed  to  be  mutually 

statistically  independent.  It  is  believed  that  assumption  3(a)  guarantees 
mutual  statistical  independence.  The  proof  together  with  simulation  results 
supporting  this  statement  is  given  in  Appendix  I.  In  Appendix  II  the 

k 

quantities  P.  are  expressed  in  terms  of  the  quantities  p0  and  it  is  proved 

k 

that,  given  the  P.  ,  the  probabilities  p0  are  zeroes  of  the  polynomial, 

J  > 

f  (A) ,  of  degree  K  in  A  given  by 
K. 


£k(X) 


AK  +  a 


K-1,K 


A^1  * 


+  ai,KX  +  a0,K 


where  the  coefficients  can  be  expressed  as 


(14) 


a 


r,K 


j=0 


(15) 


and  (.)  denotes  the  combinatorial  symbol. 

The  following  example  illustrates  this  result. 

Example  1  (K  =  3)  Referring  to  Appendix  II 

Po  >3  =  Po  po  po 

Pi  ,3  =  (1  -  Po)  Po  Po  +  (1  -  Po)  Po  po  +  (1  -  Po)  Po  Po 
P2J3  =  (1  “Po)  (1  -po)  Po  +  (1  -po)  (1  -po)  Po  +  (1  -  Po)  (1  -  Po)  Po 

Simple  algebraic  manipulation  yields 

Po  Po  Po  =  Po  ,3 

Po  Po  +  Po  Po  +  Po  Po  =  Pi  ,3  +  3P0  ,3 
Po  +  Po  +  Po  =  P2  ,3  +  2Pi  ,3  +  3P0,3 


ERL-0132-TR 


6 


Recall  that  the  coefficients  of  a  polynomial  with  zeroes  po ,  po  and  po  can  be 
computed  by  multiplying  out  the  expression 

(X  -  Po)  (X  -  Po)  (X  -  po) 

=  X3  -  (pj  +  pg  +  po)  X2  +  (p£  p l  +  pi  pi  +  po  po)  X  -  pi  pi  pi 
This  verifies  the  dense  environment  theorem  of  Appendix  II  for  K  =  3,  viz.. 


ao  ,3  =  -Po  ,3 

ai  ,3  =  Pi  ,3  +  3Po  ,3 

3-2,3  =  “  C?2  >3  +  2P,  ,3  +  3Po,a) 


k  k 

Each  emitter  has  a  unique  z  and,  once  r  is  specified,  a  unique  p0  from 

equation  (11) .  The  significance  of  the  above  result  is  that  in  a  dense 

environment  each  possible  combination  of  emitters  which  can  be  present  will 

have  a  unique  set  of  probabilities  P  ,  j  =  0,1,..., K.  The  dense  environ- 

J  >  k 

ment  of  emitters  is  said  to  be  identifiable  by  observing  the  sequence  x  . 

From  reference  4,  p392.  Ex  12.23  the  maximum  likelihood  estimate  of 
P.  ,  denoted  P.  ,  is  given  by 

J  j  ,K 


(16) 


where  P  is  given  in  equations  (6)  to  (8) 


Using  the  invariance  property 


of  maximum  likelihood  estimates  (reference  5,  Ch  VII,  Section  2.2)  the 

k  ~k 

maximum  likelihood  estimates  of  the  p0 ,  denoted  po ,  are  zeroes  of  the  poly¬ 
nomial  f  (X)  of  equation  (14)  with  coefficients 
K. 


a 


r,K 


(-D 


K-r 


r 


j=0 


(17) 


Again  using  the  invariance  property,  the  maximum  likelihood  estimate  of  z  , 
denoted  z  ,  is  given  from  equation  (12)  by 


Hk 

z 


r/[(Wk  +  1)  -  pJl 


(18) 


For  discussion  of  the  properties  of  maximum  likelihood  estimates  see 

reference  6,  Ch  4.  In  equation  (18)  Wk  is  assumed  to  be  known.  This  can 
usually  be  arranged  if  a  rough  idea  of  the  PRI  range  to  be  encountered  is 

^  1c  1c 

known.  For  example,  if  r  <  z  then  W  =0. 
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5.  EMITTERS  WITH  VARIABLE  PRI 

5.1  Single  emitters 

With  variable  PRI  emitters  the  possible  number  of  states  for  can  be 

greater  than  two.  The  number  of  states  which  can  occur  depends  on  the  type 
of  emitter  and  the  gate  time,  r.  The  probability  that  state  i  occurs  will 
be  denoted  p. . 

ri 

In  Tables  2  and  3  the  statistical  properties  of  the  observed  sequence,  x^, 

are  given  for  two  position  stagger  and  three  position  stagger  emitters 
respectively.  The  derivation  is  similar  to  that  for  fixed  PRI  emitters. 

The  formulae  in  Table  3  apply  to  both  the  123  and  132  permutations  of  PRI 
order. 

The  results  obtained  for  fixed  and  staggered  PRI  emitters  can  be  applied 
to  deriving  the  statistical  properties  of  the  observed  sequence,  x  ,  for 

emitters  with  PRIs,  y  ,  which  are  independently  and  identically  distributed 

samples  of  a  random  variable  taking  values  between  a  and  b  (a  <  b)  with 
probability  density  function  p(y  ). 

The  range  of  gate  times  is  restricted  to  r  <  3a.  Also,  the  width  of 
the  distribution  of  PRIs  is  limited  by  assuming  that  3a  >  2b.  This  implies 
2a  >  b  which  means  that  the  number  of  states  which  can  occur  will  be  three 
or  two  depending  on  r .  More  general  results  can  be  derived  but  the  required 
analysis  becomes  tedious  and  is  omitted  here. 

The  possible  number  of  pulses  which  can  occur  during  the  gate  interval  is 
shown  in  Table  4  for  the  range  of  t  of  interest  together  with  the  formulae 
for  the  state  probabilities,  p^.  These  probabilities  can  be  determined 

using  the  results  of  Table  1  for  r  <  2a  and  the  results  of  Table  2  for 
2a  <  t  <  3a  and  averaging  over  the  PRI  range  a  to  b.  The  results  of  Table  2 
must  be  used  for  the  higher  gate  time  range  because  pairs  of  pulse  repetition 
intervals  must  be  considered  for  this  case.  To  form  the  averages  the 
probability  density  function  y  p  (y)  for  r  <  2a  and  (y!  +  y2)  p(yi) 

P(y2)T2-1  for  2a  < r  <  3a  must  be  used,  where 


ri  =  I  p(y)y  dy 

a 


,b  b 

t2  =  l  I  p(yi)p(ya)(yi  +  y2)dyi  dy2 


a  ”  a 


(19) 


(20) 


t  2  =  2r  i 


(21) 


Equation  (21)  follows  from  the  fact  that 


f  p(yi)dyi 


=  i 


(22) 


The  use  of  these  probability  density  functions  to  form  the  averages  takes 
into  account  the  fact  that  a  shorter  PRI  occupies  fewer  gate  intervals  than 
a  longer  PRI  on  the  average.  Note  that  when  two  states  only  can  occur 
equivalent  results  are  obtained  if  it  is  assumed  that  a  fixed  PRI  emitter  is 
present  with  PRI  equal  to  the  average  PRI,  Tj ,  of  the  actual  emitter. 
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The  results  of  Table  4  are  now  specialized  to  two  commonly  occurring  types 
of  PRI  variation.  Upon  substituting  the  appropriate  function  p(y)  the  p^ 

for  uniform  PRI  emitters  given  in  Table  5  are  obtained.  The  corresponding 


result  for  sinusoidal  PRI  emitters  and  a  <  r  <  b  is 

po  =  I  ir  1  Ai  +  (r i  -  r)  (0.5  -  it  1  arsin  A^Tx”1  (23) 

Pi  =  [  -27T  1  Ai  +  (2r i  -  r)(0.5  +  tt  1  arsin  A2) 

+  r(0.5  -  7T_1  arcsin  A2)  j  t1~1  (24) 

p2  =  [  tt"1  Ai  +  (r  -  )  (0.5  +  7r-1  arcsin  A2)|r1"1  (25) 

where 

Ax  =  J~[  V2  -  (r  -  r  x  )2  }  (26) 

A2  =  (r  -  T1)/V  (27) 

t i  =  (a  +  b)/2  (28) 

V  =  (b  -  a)/ 2  (29) 


In  Table  6  the  state  probabilities  p0  and  p2  are  given  for  four  types  of 
emitter  PRI  at  r  =  rx  and  r  =  2rt  and  for  r  <  a.  These  probabilities  form 
the  basis  of  the  algorithm  used  to  decide  the  type  of  emitter  PRI  variation 
and  estimate  the  parameters  rx  ,  z,  etc  from  the  relative  frequencies  of 
occurrence  of  the  states  of  the  observed  data. 

An  algorithm  for  PRI  recognition  and  parameter  estimation  is  shown  in 
figure  2.  The  class  of  emitter  is  assumed  to  be  either  fixed,  two  position 
stagger  or  uniform  PRI .  The  algorithm  first  decides  the  type  of  PRI 
present  and  then  estimates  the  appropriate  PRI  parameters.  This  is 
achieved  by  observing  the  relative  frequencies  of  occurrence  of  the  states 
at  selected  values  of  t .  Parallel  processing  could  possibly  be  used  to 
obtain  the  required  data  at  the  different  gate  times.  Provided  that  the 
correct  PRI  type  is  recognized  the  parameter  estimates  obtained  can  be  shown 
to  be  maximum  likelihood  estimates. 

The  algorithm  can  be  generalized  for  other  PRI  variations.  Some 
difficulty  may  be  experienced  in  distinguishing  between  uniform  and 
sinusoidal  PRIs.  For  M-position  staggers  po  =  p2  =  0  when  r  =  Mr)  . 

This  fact  can  be  used  to  distinguish  between  complex  PRI  staggers. 

5.2  Dense  environment 

In  a  dense  environment  of  K  emitters  the  observation  xt  is  again  related 

to  the  observations  for  the  individual  emitters  comprising  the  environment 
k  . 
xt  via 


K 


k=l 
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The  set  of  states  which  can  occur  now  depends  on  r.  The  probability  that 

the  dense  environment  process  is  in  state  j,  P.  ,  can  be  related  to  the 

1 

quantities  p. .  A  polynomial  f(X)  defined  by 


fL(X) 


a.  +  a 


L-1,L 


X1-1  ♦ 


+  al,LX  *  a0,L 


(31) 


can  be  used  to  represent  a  given  environment  of  emitters  where  the 
coefficients  are  given  by 


a 


r,L 


c-uL-r 


r 


j=0 


(32) 


The  quantities  p.  are  available  upon  finding  the  factors  of  polynomial  f T (X) . 

1  L 

k 

If  r  is  less  than  the  minimum  a  then  L  =  K  and  the  analysis  of  Section  4.1 

k  k 

is  applicable  with  Tj  replacing  z  . 

It  is  possible  to  generalize  the  analysis  for  other  values  of  r .  If  for 
emitter  k  the  value  of  r  gives  rise  to  three  states  the  linear  factor 

k 

X  -  p0  in  the  polynomial  of  equations  (31)  and  (32)  is  replaced  by  the 

k  k  k 

quadratic  factor  X2  -  X(2p0  +  pi)  +  po  and  L  increases  by  1.  By  observing 
the  relative  frequencies  of  occurrence  of  the  states  it  is  possible  to 

estimate  the  polynomial  coefficients.  If  the  polynomial  is  factored 
k  k 

estimates  of  po  and  pj  can  be  obtained  which  allow  the  emitter  type  to  be 

k  k  k  k 

determined  and  the  parameters  t1  ,  z  ,  a  and  b  estimated. 


6.  ADAPTIVE  IDENTIFICATION 


Suppose  that  no  emitters  are  present.  It  is  required  to  identify  a  number 
of  emitters  which  will  then  start  transmitting.  It  is  assumed  that  when  the 
environment  of  emitters  changes  only  one  emitter  will  start  or  stop  transmitting 
and  the  time  between  changes  to  the  environment  is  long  compared  with  a 
specified  observation  interval.  Operating  under  these  conditions  the 
computational  requirements  for  emitter  identification  using  PRI  can  be  greatly 
reduced  and  improved  accuracy  achieved.  The  reason  for  this  is  that  the 
identification  method  of  this  paper  requires  a  given  polynomial  to  be  factored. 
Normally  this  entails  a  considerable  amount  of  computational  effort.  However, 
if  only  changes  to  the  environment  need  be  determined  the  problem  simplifies 
because,  in  splitting  a  polynomial  into  two  factors,  the  second  factor  is 
easily  found  once  the  first  factor  is  specified.  This  approach  will  be  called 
adaptive  identification. 

Suppose  that  a  number  of  emitters  are  present  and  the  polynomial  associated 
with  this  environment  is  given  by 


fLw  -  xL  *  Vi xL'1 


(33) 


+  ...  +  ai  X  +  ao 
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(In  this  section  a  is  replaced  by  a  to  simplify  the  notation) . 

IT  2  Li  r 

If  another  emitter  starts  transmitting  prior  to  the  next  observation  interval  and 
for  this  emitter  two  states  are  observed  at  the  value  of  r  selected,  then  the 
new  environment  can  be  represented  by  the  polynomial 

f2  (X)  =  XL+1  +  b.  XL  +  ...  +  b,  X  +  b0  (34) 

L+ 1  L 

It  is  easy  to  show  that  the  state  probability,  p<> ,  for  this  emitter  is  given  by 

p°  =  aL_1  -  bL  =  -b0/a0  =  (ai  ;1  -  bi)/ai 

i  =  1,2,...,  L-l  (35) 

Thus,  it  is  possible  to  estimate  the  new  emitter  parameter  with  simple  rational 

operations  if  the  prior  knowledge  about  the  environment  is  used.  The  coefficients 

a.  are  assumed  known  and  the  coefficients  b.  are  estimated  from  observations  of 
i  l 

the  relative  frequencies  of  occurrence  of  the  states. 

If  the  new  emitter  has  three  states  at  the  value  of  r  selected  the  new 

environment  is  represented  by  the  polynomial 


f2 

L+2 

CM 

=  XL  ^  +  b^+  ^  XL  ^  +  ...  +  biX  +  bo 

(36) 

The  state  probabilities 
coefficients  via 

p0  and  p!  for  this  emitter  are  related  to 

the  polynomial 

Cl 

= 

-(2p0  +  pi) 

(37) 

Co 

= 

Po 

(38) 

bo 

= 

Co  So 

(39) 

bi 

= 

Cl  So  +  Co  ai 

(40) 

b. 

i 

= 

ai  2  +  ^i  +  Co a^ }  x  —  2 1  ••>,  L  -  1 

(41) 

bL 

= 

aL-2  +  ClVl  +  Co 

(42) 

bL+l 

= 

aL-l  +  Cl 

(43) 

A  theoretical  framework  for  deriving  optimal  adaptive  parameter  estimates  is 
not  known  at  this  stage.  Further  discussion  of  adaptive  identification  is 
postponed  until  Section  8. 


11 


ERL-0132-TR 


7.  PULSE  DETECTION 

In  order  to  identify  radar  emitters  using  PRI  the  number  of  pulses  occurring 
within  the  gate  interval  must  be  determined.  If  the  pulses  do  not  overlap  it 
is  relatively  easy  to  determine  the  number  of  pulses  which  have  occurred.  In 
practice,  the  probability  that  pulses  will  overlap  is  small  and  it  may  be 
possible  to  assume  that  overlap  never  occurs.  Nevertheless,  in  this  section  the 
possibility  of  detecting  overlapping  pulses  will  be  investigated. 

Suppose  that  the  received  RF  signal,  rRp(s),  from  K  emitters  is  represented 

at  the  time  s  by 

K 

A^(s)  cos^f^s  +  (44) 

k=l 

k 

where  the  pulse  amplitude  modulation  waveform  is  denoted  by  A  (s) .  The  radio 

k  k 

frequency  and  phase  of  emitter  k  are  denoted  by  f  and  y?  respectively. 

Let  the  centre  frequency  of  the  RF  filter  be  f<> .  The  IF  signal  in  the 
receiver  can  be  expressed  as 


rRF(s) 


riFCs)  =  Aip(s)  cos[  271’fjpS  +  aiF(s)j 


(45) 


where  the  amplitude,  A  (s) ,  and  phase,  a  (s) ,  are  given  by 

1  r  1  r 


Aif(s)  = 


K  K 


k  „r,  k  r? 


k=l 


fA^(s)]2  +  A^(s)Ar(s)  cosl27r(fK  -  fr)  s  +  -  <p 

k=l  r=l 
k  ^  r 


(46) 


(iIF(s)  =  tan-1 


A^(s)  sin^2ir(f^  -  f0)s  + 


k=l 


^ ^  Ak(s)  cos[27r(f^  -  f0)s  + 


k=l 


(47) 


The  derivation  of  equations  (45)  to  (47)  is  discussed  in  more  detail  in 
reference  7.  In  what  follows  the  radar  pulses  received  from  each  emitter  are 
assumed  to  be  approximately  rectangular  i,n  shape. 

If  f  ^  f  for  k  r  then  the  time  average  of  AT„(s)  will  increase  sharply 

whenever  a  new  pulse  occurs  no  matter  if  other  pulses  are  already  present. 

If  f^  =  fr  for  some  k  #  r  then  except  under  freak  conditions  the  arrival  of  a 
new  pulse  will  be  accompanied  by  a  sharp  change  in  the  average  phase  or 
amplitude  of  the  IF.  Similarly,  when  a  particular  pulse  finishes  the  IF  phase 
and  amplitude  will  change  sharply. 
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From  the  preceding  analysis  it  would  appear  that  it  is  possible  in  principle 
to  detect  overlapping  pulses.  The  practical  implementation  of  a  suitable 
detector  could  present  problems.  The  use  of  Bragg  cell  technology  can  avoid  all 
these  problems  with  detecting  overlapping  pulses.  At  present  these  devices 
are  relatively  expensive  and  limited  in  dynamic  range.  Two  alternate  approaches 
are  considered  next. 

It  may  be  possible  to  develop  a  PRI  measurement  system  where,  after  a  pulse 
has  been  detected,  all  further  signal  activity  is  ignored  for  a  specified  dead 
time.  The  dead  time  would  be  comparable  to  the  longest  expected  pulse  width, 
say  10  Ms.  The  statistical  analysis  of  this  paper  would  need  to  be  modified 
to  take  this  dead  time  into  account.  The  equations  would  be  a  little  more 
complicated  but  the  problem  of  detecting  overlapping  pulses  could  then  be 
ignored. 

Another  approach  is  to  use  bandpass  filters  and  logic  gates.  Consider  the 
following  simple  illustrative  example.  Suppose  that  signal  activity  at  two 
different  RF  values  is  detected.  Two  bandpass  filters  with  pulse  detectors 
connected  to  an  exclusive  OR  gate  as  shown  in  figure  3  could  be  used  to  facilitate 
PRI  analysis.  Provided  that  exactly  coincident  pulses  with  identical  pulse 
width  do  not  occur  (a  low  probability  event)  overlapping  pulses  will  enter  the 
pulse  counter  as  two  well  defined  pulses.  The  PRI  parameter  measurement  method 
proposed  in  this  paper  could  then  be  applied  using  the  data  obtained  with  this 
simple  device. 

If  missing  pulses  occur  (or  pulses  in  a  train  are  not  detected)  then  the  PRI 
measurements  will  of  course  be  degraded.  Precautions  would  need  to  be  taken  to 
avoid  errors  caused  by  missing  pulses.  Detailed  discussions  of  this  topic  are 
left  for  later  work. 


8.  EXPERIMENTAL  RESULTS 

In  this  section  computer  simulations  and  laboratory  work  which  have  been 
performed  to  test  the  emitter  PRI  measurement  method  of  this  paper  are  described. 

8.1  Fixed  PRI  emitters 

The  results  presented  in  figures  4  and  5  demonstrate  the  accuracy  which 

can  be  achieved  and  show  the  importance  of  the  choice  of  r .  If  the  maximum 

likelihood  estimates  are  accurate  then  the  zeroes  of  the  polynomial  of 

£quation  (14)  will  be  real.  For  short  observation  intervals  the  estimates, 

P.  ,  may  be  inaccurate  resulting  in  complex  zeroes  for  the  polynomial. 

J  j  k 

In  this  case  the  PRI  estimation  method  can  not  be  used  as  equation  (18)  is 
no  longer  valid.  To  investigate  this  problem  simulations  were  carried  out 
with  K  =  2,  z1  =  103  Ms  and  z2  =  223  Ms  (nominally).  In  figure  4  the 
percentage  of  runs  for  which  complex  zeroes  occur,  CR,  is  plotted  versus  r 

for  three  observation  intervals  of  duration  0.01,  0.1  and  1.0  s.  The 
percentage  of  complex  zeroes  decreases  for  longer  observation  times  and  is 
considerably  reduced  for  certain  values  of  r.  In  figure  5  a  performance 
index  P^  to  be  defined  next  is  plotted  versus  r  for  the  same  three 

observation  intervals.  The  performance  index  is  given  by 

Pj  =  2(1  -  CR/100)  (E1  +  E2)-1  (48) 

k  -4c 

where  E  is  the  r.m.s.  (root  mean  square)  percentage  error  in  z  over  a 

k  k 

number  of  runs  with  d0  sampled  from  a  uniform  distribution  on  0  to  z  . 

The  number  of  runs  is  104 ,  103  and  102  for  the  0.01,  0.1  and  1.0  s  observ¬ 
ation  intervals  respectively.  The  greater  the  value  of  this  index  the 
nearer  the  PRI  estimates  are  expected  to  be  to  the  true  PRIs.  The 
estimates  may  be  deemed  to  be  acceptable  if  P  >  1.  Notice  that  certain 
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values  of  r  yield  more  approximate  PRI  estimates. 

In  figure  6  the  discriminant  of  the  quadratic  function  f2  (A) , 

(ai  ,1  -  4  ao,j),  is  plotted  versus  r.  Referring  to  the  results  of  figures  4 
and  5  it  can  be  seen  that  improved  performance  is  achieved  for  values  of  r 
for  which  the  discriminant  of  the  quadratic  function  is  greater. 

It  may  be  possible  to  select  r  to  obtain  optimal  discrimination  between 
emitters  if  the  scenario  to  be  encountered  is  known.  However,  in  many  cases 
a  fixed  value  of  r,  say  100  Ms,  will  be  necessary  due  to  practical  limitations. 
This  case  is  examined  next. 

The  results  obtained  in  the  figure  7  show  the  accuracy  which  can  be 
obtained  with  a  fixed  value  of  t.  The  graphs  in  figures  7(a)  and  7(b)  are 
of  the  percentage  of  runs  for  which  complex  zeroes  occur  (K  =  2)  and  r.m.s. 
percentage  error  (K  =  1  and  K  =  2)  respectively  versus  observation  time. 

The  value  of  r  is  100  Ms.  The  number  of  runs  is  1000  and  for  each  run  the 

PRIs  were  selected  randomly  from  the  range  200  to  5000  Ms.  When  K  =  2  the 

zeroes  of  a  quadratic  are  determined  to  obtain  PRI  estimates  for  the  two 

emitters  present.  As  can  be  seen  the  results  are  not  very  good.  A  large 

error  results  if  only  two  emitters  are  present  and  complex  zeroes  occur  quite 
frequently,  which  renders  the  method  almost  useless.  These  problems  could 
be  overcome  by  altering  r  from  its  100  Ms  value.  Perhaps  a  better  way  is 
to  use  the  adaptive  approach  of  Section  6. 

Recall  that  with  the  adaptive  approach  all  the  factors  of  the  polynomial 
representing  the  environment  need  not  be  found  as  the  change  to  the 
environment  can  be  deduced  simply  by  comparing  polynomial  coefficients. 

Prior  knowledge  of  the  environment  allows  the  coefficients,  a^,  to  be 

precisely  calculated.  The  coefficients,  b^,  of  the  polynomial  associated 

with  the  new  environment  are  estimated  using  the  relative  frequencies  of 
occurrence  of  states  of  the  observed  data.  Two  estimates  of  PRI  given  by 
equations  (50)  and  (51)  are  considered 


Z  =  T /[  (W  +  1)  -  Poj 


(49) 


po  =  (L  +  1) 


-i 


L-l 

*l-i  -  bl  -  Wa» +  y  (ai-i  ■  ®i)/ai 

i=l 


(50) 


po  =  -  bL 


(51) 


where  L  is  the  number  of  fixed  PRI  emitters  present  prior  to  the  transmission 
of  the  new  emitter. 

The  results  for  the  estimates  of  equations  (50)  and  (51)  are  shown  in 
figures  8  and  9  respectively.  The  number  of  runs,  the  value  of  r  and  the 
PRI  ranges  considered  are  given  on  the  figures.  A  curve  is  labelled  L 
emitters  present  if  (L  -  1)  emitters  with  precisely  known  PRI  comprise  the 

environment  prior  to  the  transmission  of  the  L*"*1  emitter,  whose  PRI  is  to  be 
estimated.  The  r.m.s.  %  error  in  the  resulting  estimate  is  plotted.  It 
is  seen  that  the  results  with  the  estimate  of  equation  (51)  are  superior  and 
for  up  to  five  emitters  present  and  one  second  observation  times  the  r.m.s. 
percentage  error  in  the  estimates  obtained  is  less  than  0.4%. 
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8.2  Variable  PRI  emitters 

The  generalisation  of  the  adaptive  identification  approach  to  variable 
PRI  emitters  is  discussed  now. 

If  r  is  such  that  for  the  new  variable  PRI  emitter  only  two  states  are 
observed  then  the  approach  of  Section  8.1  can  be  used  to  estimate  the  average 
PRI,  t j  ,  with  t'i  replacing  z  in  equation  (49).  If  r  is  such  that  for  the 
new  variable  PRI  emitter  three  states  are  observed  p0  and  pi  are  estimated 
by  pb  and  pi  given  as 


Po 

=  Co 

(52) 

Pi  = 

-Ci  -  2 Co 

(53) 

Cl  = 

®L+1  “  aL-l 

(54) 

Co  =  bL  "  aL-2  '  Cl  Vl  (55) 

where  the  degree  of  the  polynomial  representing  the  initial  environment  is  L. 
These  estimates  are  selected  because  it  appears  that  estimates  are  best 
calculated  using  the  highest  order  coefficients  in  the  polynomial.  The 
estimates  p^  and  Tj  can  be  used  in  the  algorithm  of  Section  5.1  to  decide 

the  type  of  variable  PRI  for  the  new  emitter  and  estimate  the  parameters 

t i  ,  z  ,  a  and  b  of  the  new  emitter, 
m 

Simulation  results  for  adaptive  identification  applied  to  variable  PRI 

emitters  are  given  in  figures  10  and  11.  Background  environments  of  0,  1, 

2  and  3  emitters  are  considered.  The  emitters  for  the  background 

environment  are  selected  from  a  614  Ms  fixed,  a  401/491  Ms  2  position 

stagger  and  a  950  to  1050  Ms  uniform.  The  new  emitter  can  be  either  a 

761  Ms  fixed,  a  721/887  Ms  2  position  stagger  or  a  813  to  997  Ms  uniform. 

For  a  given  background  environment  the  type  of  PRI  of  the  new  emitter  is 

determined  and  then  the  parameters  ,  z  a  and  b  are  estimated.  In 

m 

figure  10  the  probability  of  error  in  determining  the  emitter  type  is 
plotted  versus  observation  time  for  each  of  the  three  new  emitters  considered. 
The  algorithm  of  figure  2  is  used  with  the  decision  thresholds  er  and  e2 
shown  on  figure  10.  The  r.m.s.  percentage  error  in  the  estimates  of  r i ,  z^ 

a  and  b  are  plotted  versus  observation  time  in  figure  11  assuming  that  the 

new  emitter  is  correctly  classified  in  each  case.  These  results  show  that 

the  type  of  emitter  PRI  can  be  recognized  and  parameter  estimates  subsequently 
extracted  using  the  techniques  proposed  in  this  paper. 

The  implementation  of  these  concepts  using  both  analogue  and  digital 
techniques  is  being  investigated.  The  digital  techniques  show  the  greatest 
promise.  The  system  under  development  consists  of  a  pulse  detector,  a  4-bit 
counter  and  an  Intel  8080  microprocessor.  The  environment  of  emitters  for 
testing  the  identification  system  is  provided  by  a  microwave  laboratory 
simulator.  The  4-bit  counter  is  used  to  count  the  number  of  pulses  which 
occur  during  the  gate  interval.  The  gate  time  is  controlled  by  the  micro¬ 

processor.  More  details  will  be  given  in  a  later  report. 
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9.  DISCUSSION 


9.1  Introduction 

Radar  emitter  recognition  in  a  rapidly  changing  dense  environment  of 
emitters  with  parameter  agility  is  a  challenging  problem. 

"One  signal  processing  area  that  needs  immediate  improvement  is  the 
handling  of  complex  time  domain  modulations.  IFMs  and  other  wide  open 
receivers  can  handle  RF  agility,  but  time  domain  agility  remains  a  serious 
problem  for  all  receiver  types.  It  is  a  very  difficult  task  to  correctly 
process  a  signal  having  multiposition  stagger  or  PRI  shifts  in  the  presence 
of  other  emitters  without  generating  false  emitters  or  incorrect  IDs. 

Angle  of  arrival  is  a  powerful  tool  in  the  solution  of  this  problem (ref. 8)". 

This  quotation  testifies  to  the  importance  of  the  signal  processing 
problem  addressed  in  this  paper.  The  suggestion  of  directional  antennas 
may  not  always  provide  a  solution.  It  may  not  always  be  possible  to  thin 
the  environment  sufficiently  using  direction  of  arrival.  The  limitations 
of  cost,  space,  time  and  higher  priority  tasks  such  as  direction  finding 
and  directive  countermeasures  have  to  be  taken  into  account  also. 

Several  recognition  methods  using  PRI  measurement  are  adopted  in  practice 
and  this  paper  puts  forward  yet  another  approach.  Each  method  relies  on 
the  validity  in  the  real  world  of  the  particular  assumptions  involved. 

Other  methods  of  radar  emitter  PRI  sorting  are  described  in  this  section 
and  then  some  features  of  the  method  proposed  here  are  discussed. 

9.2  Other  methods 

Current  methods  of  PRI  analysis  are  mainly  applicable  to  emitters  with 
fixed  PRI.  The  pulse  train  sorter  passes  the  received  pulse  train  and  a 
delayed  version  through  a  coincidence  circuit.  If  the  delay  is  equal  to 

the  emitter  PRI  an  output  is  obtained.  Another  method  uses  the  times  of 
arrival  of  the  incoming  pulses.  The  time  interval  between  the  next  two 
unsorted  pulses  is  measured  and  then  if  pulses  are  found  at  2,  3,  4  ...,  r 
times  this  interval  from  the  first  pulse  an  emitter  with  this  PRI  is  deemed 
to  be  present.  Depending  on  the  search  method  employed  and  the  frequency 
with  which  missing  pulses  occur  both  these  methods  are  prone  to  respond  to 
subharmonic  pulse  repetition  frequencies. 

Calculating  the  autocorrelation  function  of  the  incoming  pulse  train 
is  also  suggested  as  a  possible  means  of  PRI  sorting.  For  a  single  emitter 
with  fixed  PRI  the  autocorrelation  function  bears  a  striking  resemblance  to 
the  original  pulse  train.  Calculating  the  autocorrelation  function 
digitally  would  involve  a  considerable  amount  of  computational  effort,  but 
optical  signal  processing  techniques  may  allow  the  autocorrelation  function 
to  be  calculated  almost  instantaneously.  Unfortunately,  a  systematic 
method  of  analysing  the  autocorrelation  function  for  a  dense  environment  of 
variable  PRI  emitters  does  not  appear  to  be  available  and  handling  the  data 
rate  involved  would  not  be  a  trivial  task. 

A  common  problem  with  these  standard  methods  is  that,  if  implemented  on 
a  microprocessor,  the  processor  rapidly  becomes  overloaded  in  a  dense 
environment.  Even  if  elaborate  buffering  is  provided  it  may  be  necessary 
to  ignore  some  received  data. 

9.3  Blanking 

Blanking  is  a  technique  for  thinning  the  environment.  Time  domain 
blanking  can  be  used  with  all  methods  of  emitter  identification.  After  an 
emitter  is  identified  further  pulses  which  are  deemed  to  have  come  from  this 
emitter  are  blanked  or  ignored  in  the  analysis  of  future  received  data. 

This  approach  is  quite  successful  with  emitters  with  fixed  PRI  in  a  dense 
environment  where  the  emitters  appear  one  at  a  time  and  at  a  rate  which  does 
not  saturate  the  system. 
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Time  domain  blanking  could  well  be  applied  to  boost  the  number  of  emitters 
which  can  be  successfully  identified  in  a  dense  environment  using  the  PRI 
measurement  method  proposed  in  this  paper. 

9.4  Scanning  antennas 

A  requirement  of  the  PRI  measurement  method  proposed  in  this  paper  is 
that  the  environment  does  not  change  rapidly.  Sufficient  data  must  be 
collected  while  the  environment  is  static  to  obtain  accurate  PRI  estimates. 

If  emitters  with  scanning  antennas  are  present  and  only  main  lobe  pulses 
can  be  detected  this  requirement  may  not  be  met.  Possible  solutions  to 
this  problem  are  discussed  in  what  follows. 

To  increase  the  number  of  observed  pulses  in  dealing  with  search  radars 
the  threshold  of  the  pulse  detection  circuit  could  be  lowered  to  detect  pulses 
in  the  antenna  side  lobes.  The  drawback  with  this  is  that  it  may  preclude 
emitter  sorting  using  other  parameters  such  as  direction  of  arrival  and 
radio  frequency.  Also  the  scan  information  for  the  emitter  concerned  will 
be  lost  and  in  a  noisy  environment  the  resulting  false  alarm  rate  could  be 
too  high.  To  achieve  the  required  receiver  sensitivity  a  narrow  bandwidth 
would  be  necessary.  This  would  no  doubt  reduce  the  probability  of  intercept. 

A  better  approach  is  suggested  now.  Suppose  that  data  are  collected 
using  a  series  of  consecutive  observation  intervals  each  of  which  is  shorter 
than  the  shortest  main  lobe  pulse  burst  to  be  encountered.  Accurate  PRI 
estimates  could  not  be  obtained  using  data  from  any  individual  observation 
interval.  However,  it  should  be  possible  to  determine  when  the  main  lobe 
of  a  particular  emitter  occurs  by  correlating  the  data  from  all  the 
consecutive  observation  intervals.  By  averaging  the  data  recorded  during 
successive  main  lobe  bursts  it  is  conjectured  that  accurate  PRI  parameter 
measurements  can  be  obtained  for  emitters  with  scanning  antennas.  The  data 
from  observation  intervals  in  between  main  lobe  pulse  bursts  will  be  useful 
too  if  nonscanning  radar  emitters  are  present. 

9.5  Very  short  PRI 

In  practice  the  gate  time  can  not  be  made  arbitrarily  short  because  a 
certain  amount  of  processing  must  be  done  during  this  time.  PRI  of  the 
order  of  a  few  microseconds  can  occur  for  some  emitters.  This  problem  can 
be  dealt  with  by  treating  such  emitters  as  a  special  case.  They  will 
produce  a  recognizable  factor  in  the  polynomial  representing  the  environment 
for  any  given  value  of  r . 

9.6  Recognition  of  known  emitters 

Often  it  is  necessary  to  be  able  to  quickly  recognize  that  a  particular 
emitter  has  started  transmitting.  This  is  the  case  with  certain  threat 
emitters.  The  polynomial  factor  associated  with  such  an  amitter  could  be 
stored  in  a  library  and  the  adaptive  method  of  this  paper  used  to  detect 
the  presence  of  this  emitter.  This  approach  would  be  advantageous  for 
threats  with  PRI  agility. 

9.7  Sorting  using  pulse  amplitude 

The  adaptive  technique  for  PRI  measurement  presented  in  this  paper  examines 
changes  in  the  environment.  It  is  possible  to  artificially  obtain  the 
effect  of  an  emitter  starting  or  stopping  transmitting  using  pulse  amplitude 
selection.  Sorting  using  pulse  amplitude  in  conjunction  with  the  adaptive 
PRI  measurement  technique  would  appear  to  be  a  promising  method  for  emitter 
identification. 

9.8  Summary 

The  PRI  measurement  method  proposed  here  can  be  viewed  as  a  generalization 
of  the  familiar  pulse  counter  to  deal  with  variable  PRI  in  a  dense  environ¬ 
ment.  The  method  is  adaptive  in  that  it  analyses  changes  to  the  environment 
to  obtain  parameter  measurements.  One  advantage  of  the  method  is  its 
simplicity  which  makes  it  ideal  for  implementation  on  a  microprocessor. 
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It  does  not  involve  time  consuming  parameter  searches  even  to  deal  with 
variable  PRI  emitters.  It  allows  emitters  to  be  identified  with  a  high 
probability  of  intercept.  Further  work  is  in  progress  to  implement  the 
concept  using  microprocessor  technology. 


10.  CONCLUSION 

The  technique  described  in  this  paper  should  have  applications  to  surveillance 
problems  in  electronic  warfare.  Although  alternate  techniques  exist,  the 
method  of  this  paper  allows  parameter  agile  radar  emitters  to  be  identified  in  a 
dense  environment  comparatively  cheaply  and  simply. 
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NOTATION 

For  the  following  symbols  if  a  superscript,  k,  is  added  this  denotes  that  the 

symbol  is  associated  with  emitter  k.  Also,  if  ~  is  placed  over  a  symbol  this 

indicates  an  estimate  of  the  quantity  denoted  by  the  symbol. 

a  parameter  of  a  variable  PRI  emitter  specifying  PRI 

distribution 

a  .  coefficient  r  in  f.(X) 

r,j  3V  J 

a  coefficient  r  in  f1.  (X) 

r  3 v 

b  parameter  of  variable  PRI  emitter  specifying  PRI  distribution 

b  coefficient  r  in  f2.  (X) 

r  3 

Co  coefficient  of  quadratic  factor  defined  in  equation  (37) 

Ci  coefficient  of  quadratic  factor  defined  in  equation  (38) 


fj(*) 

f-  (X) 
f*  00 

fk 

fIF 


Pi 

P(yn) 

riF(s) 


z 

z 

m 


polynomial  X-1  +  a.  ,  .  X-^  *+...+  a,  .  X  +  a^  . 

3-1.3  1.3  0,3 

polynomial  of  degree  j  representing  initial  environment 

polynomial  of  degree  j  representing  environment  after  new 
emitter  starts  transmitting 

radio  frequency  of  kcn  emitter 

IF  frequency 

denotes  emitter  number 

number  of  pulses  observed  in  gate  interval  t 

probability  that  xt  =  i 

probability  density  function  of  PRI 

IF  signal  at  time  s 

RF  signal  at  time  s 

time  variable 

time  point  t 

subscript  denoting  discrete  time  points 
observation  derived  from  n^  for  gate  interval  t 
observation  at  time  point  t  when  emitter  k  only  is  present 
PRI  sample 

PRI  of  fixed  PRI  emitter 
staggered  PRI  parameter 
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At 

A2 

Ak(s) 


aif0) 

CR 

Ek 

K 


M 


N 

q 

N  . 
mm 


N 

max 


P 

q 

Q 

T 

V 

W 

aIF(s) 

8(0 

A 


r 

Tl 

r2 


|Yn,n=l,2,...,N] 


defined  in  equation  (26) 
defined  in  equation  (27) 

amplitude  of  component  of  RF  signal  due  to  emitter  k  at 
time  s 

amplitude  of  IF  signal  at  time  s 

percentage  of  simulations  for  which  complex  zeroes  occur 

root  mean  square  percentage  error  in  ^ 

number  of  emitters  present 

number  of  positions  for  PRI  stagger 

defined  in  equation  (7) 

defined  in  equation  (8) 

minimum  number  of  pulses  detected  in  any  gate  interval 
during  observation  interval 

maximum  number  of  pulses  detected  in  any  gate  interval 
during  observation  interval 

performance  index  of  equation  (48) 

probability  of  state  j  occurring  in  dense  environment 
of  K  emitters 

relative  frequency  of  occurrence  of  state  q 

number  of  observed  states  is  Q  +  1 

number  of  discrete  time  points  (or  gate  intervals) 

half  width  of  PRI  distribution 

non-negative  integer  specifying  gate  time  limits 

phase  of  IF  signal 

Dirac  delta  function 

polynomial  variable 

time  measured  from  start  of  gate  interval  to  where  a 
specified  pulse  occurs 

gate  time 

average  PRI 

twice  average  PRI 

phase  of  emitter  k 

combinatorial  symbol  r  !/f  (r-j)  !j  !} 

stochastic  process  representing  PRI 
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APPENDIX  I 

VERIFICATION  OF  STATISTICAL  MODEL 

1.1  (Section  4.1)  0  may  be  considered  to  be  a  sample  value  of  a  random  variable 
uniformly  distributed  from  0  to  z 

Assume  that  the  gate  interval  partition  and  the  radar  pulse  train  are 
synchronized  to  a  clock  of  period,  7,  which  can  be  made  arbitrarily  short. 
Thus,  r  =  n7  and  z  =  m7  where  n  and  m  are  positive  integers.  Assumption  3(b) 
then  becomes  equivalent  to  assuming  that  m  and  n  are  relative  prime. 

Also,  0  =  j  7  where  1  ^  j  <  m.  It  is  required  to  prove  that  j  is 

uniformly  distributed  on  the  set  l,2,3,...,m  .  Note  that 


3q+r  =  3q  +  (W  +  l)m  -  n  ]  (mod  m) 


where  W  is  defined  in  Table  1.  The  sequence,  [  j  j ,  is  periodic  with  period 

m.  There  are  mn  clock  periods  between  the  pulses  associated  with  j  and 
jq+m  and  each  value  of  occurs  once  during  a  period.  If  this  were  not 

true  then  for  some  q  and  i,  1  <  i  <  m,  =  Jq+i-  This  would  contradict 

the  assumption  that  m  and  n  are  relatively  prime. 

This  result  can  be  generalized  further  to  the  case  where  j  can  be  real, 

that  is  jt  e  (0,m] .  This  corresponds  to  the  case  where  the  radar  pulse 

train  is  not  exactly  synchronized  with  the  clock.  Suppose  that  for  some  t 


0t  =  (it  -  A)t 

for  some  A  e  [0,1)  and  i^  e  [  1,2, . . . ,m] .  Then  0  can  be  expressed  in  this 

form  for  all  t.  The  sequence,  [i  is  periodic  as  before  and  jt  is 

uniformly  distributed  on  [ 1  -  A,  2  -  A,  . . . ,  m  -  A] . 

For  practical  purposes  the  gate  interval  partition  and  the  radar  pulse 
train  are  unsynchronized  if,  for  some  choice  of  A  and  7,  relatively  prime  m 
and  n  can  be  found  (with  m  sufficiently  large)  which  allow  an  accurate  model 
to  be  obtained.  If  these  conditions  are  met  then  the  distribution  of  0 

will  closely  approximate  a  uniform  distribution  on  0  to  z,  independent  of  the 
initial  phase  of  the  radar  pulse  train  relative  to  the  gate  interval 
partition. 

k 

1.2  (Section  4.2)  The  sequence  of  states  x  associated  with  the  individual 

emitters  are  mutually  statistically  independent 

Assume  that  the  gate  interval  partition  and  the  radar  pulse  trains  are 
synchronized  to  a  clock  of  period,  7,  which  can  be  made  arbitrarily  short, 
k  k  k 

Thus  T  =  ny  and  z  =  m  7  where  n  and  m  are  positive  integers.  Assumptions 

j/- 

3(a)  and  (b)  then  become  equivalent  to  assuming  that  integers  n,  m1  ,  ...,  m  , 

k  k  k  k 

are  relatively  prime.  Note  that  0  =  j  7,  where  1  <  j  <m  .  Using  a 

L  u  t 

~  K.  K 

similar  argument  to  that  in  (1.1)  above  the  vector  sequence,  J*  =  [ 

K  K 

is  periodic  with  period,  m  ...  m  ,  and  each  possible  value  of  occurs 
once  for  this  number  of  successive  gate  intervals.  Thus,  the  sequence. 
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K  (  K? 

. ..,  ,  can  be  regarded  as  uniformly  distributed  on 

i  K  K  k 

(0,  z1]  x  ...  x(0,  z  ] .  Since  the  value  of  6  determines  the  value  of  x 

t  V 

this  is  tantamount  to  having  mutually  statistically  independent  sequences  of 

+  k 
states,  x  . 

1.3  Model  validation  via  simulations  for  fixed  PRI  emitters 

The  results  in  figures  12  and  13  are  submitted  as  evidence  that  the  model 
introduced  faithfully  represents  the  observed  data.  The  radar  emitter 
pulse  trains  and  the  observation  procedure  are  simulated  in  double  precision 

k 

FORTRAN  on  a  32  bit  machine  so  that  a  number  of  rims  with  varying  z  ,  r  and 

k  k 

60  can  be  performed.  For  each  run  z  and  r  are  samples  from  uniform 

distributions  with  specified  limits  and  &0  is  a  sample  from  a  uniform 

jc 

distribution  on  the  interval  0  to  z  .  Consequently,  the  probability  that 
assumptions  3(a)  and  3(b)  are  invalid  is  vanishingly  small.  The  error  for 
run  R,  E_,  is  recorded  as 

K 

K 

er  *  +  ]T 

j=0 

and  then  the  root  mean  square  percentage  error,  E,  computed  via 


For  the  simulations  carried  out  N^  =  1000.  In  figures  12  and  13  E  is 

plotted  versus  the  duration  of  the  observation  interval.  In  both  figures  r 

v 

lies  between  90  and  100  Ms.  In  figure  12  z  lies  between  103  and  303  ms 

k 

and  curves  are  plotted  for  K  =  2,  3  and  4.  In  figure  13  K  =  2  and  z  lies 

v 

between  103  and  303  Ms  for  curve  C  and  z  lies  between  103  and  1103  Ms  for 
curve  D.  The  fact  that  the  curves  in  figures  12  and  13  all  decrease  to 
fairly  small  values  indicates  that  the  model  introduced  is  valid.  Notice 
that  convergence  is  slower  when  K  is  larger  and  when  the  range  of  allowable 
PRIs  is  greater. 

1.4  Model  validation  via  simulations  for  two  position  stagger  PRI  emitters 

The  simulations  described  next  justify  the  use  of  the  proposed  models  for 
two  positions  stagger  PRI  emitters.  The  radar  emitter  pulse  trains  and  the 
observation  procedure  are  simulated  in  double  precision  FORTRAN  on  a  32  bit 

k  k 

machine  so  that  a  number  of  runs  with  varying  z  ,  r  and  do  can  be  performed, 
k 

For  each  run  z^  and  r  are  samples  from  uniform  distributions  with  specified 

k  k 

limits  and  60  is  a  sample  from  a  uniform  distribution  on  either  0  to  z  or 

k  k 

0  to  zi  +  z2  as  appropriate.  Consequently  the  probability  that  assumptions 
3(a)  and  (b)  are  invalid  is  vanishingly  small.  The  root  mean  square 
percentage  error,  E,  defined  above  is  computed  for  N^  =  1000.  The  curves 
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shown  in  figure  14  are  for  one  fixed  and  one  staggered  PRI  emitter  present 

and  for  one  fixed  and  two  staggered  PRI  emitters  present.  In  figure  14  E 

is  plotted  versus  the  length  of  the  observation  interval.  For  the  emitter 

with  fixed  PRI  z1  lies  between  103  and  303  Ms.  For  the  two  position  stagger 

k  k  lc 

PRI  emitters  Zi  +  z2  lies  between  203  and  603  Ms  and  Zj  and  r  are  selected 

with  the  following  constraints 

50  <  zi'  <  [  z!'  +  z^j  •  /2 

c  J  irmv 


max 


<  T  <  lz£} 


min 


r  <  z1 

The  fact  that  the  curves  both  decrease  to  small  values  indicates  that  the 
models  introduced  are  valid.  Convergence  is  slower  when  more  emitters  are 
present. 
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APPENDIX  II 

DENSE  ENVIRONMENT  THEOREM 


This  appendix  contains  the  proof  of  the  result  of  Section  4.2  that,  given  the 

k 

P.  ,  the  probabilities  p0  are  zeroes  of  the  polynomial,  f  (X ) ,  of  degree  K  in 

j  ,R  R 

A  given  by 


fK(X)  »  XK  *  Vl,KXK'1  *  •••  *  ai,KX  *  a0,K 


where  the  coefficients  can  be  expressed  as 


r,K 


(-1) 


K-r 


K-j 
j,K  \r- j 


j=0 


and  (.)  denotes  the  combinatorial  symbol. 


k  . 


To  simplify  notation  in  this  appendix  po  is  replaced  by  p. 


k  „ 


Using  the  statistical  independence  of  the  individual  processes  xt  discussed 
in  Appendix  I  the  expression  for  P^  ^  is 


0SK 


=  Pi  P2 


...  PK 


p 

j,K 


(1  -  p.  )(1  -  p.  )...(!  -  p.  )  p.  -  p.  2  Pi 

X1  X2  Xj  Xj  1  Xj  2  XK 


ii  ,U  ,  •  •  -,1^  =  1 


ii  ^12  ^  *  <  ij 

{  i.  1,  i .  "**  2,  .  . . ,  i._j  —  [  1*  2,  ...,  K]  —  [  ij  ,  i2  ,  .  .  • ,  i .] 

DIR  J 


,j  -  1 K . 


Another  useful  relation  is  for  1  <  j  <  K 


Pj  ,K+1  Pj,K  PK+1  +  Pj-1,K  (1  ~  Pk+P 


The  theorem  is  proved  by  means  of  the  following  three  lemmas. 
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Lemma  1:  Let  pi  ,  p2  ,  . . . ,  p„  be  zeroes  of  the  polynomial  of  degree  K  in  X 
“  K 


f„(x)  =  xK  +  av  „  „  XK_1  +  ...  +  a,  „  X  +  a„ 


V 


K-1,K 


1,K 


0,K 


(II.  1) 


That  is,  for  r  =  1,  2,  . ..,  K-l 


K 


\,K  -  <-« 


K-r 


P*  P,  P, 

ll  12  1 


ii  >  i 2  » •  •  •  1 


K-r 


ii  <  i2  <  •  .  .  <  i 


K-r 


(II. 2) 


and 


a0,K  =  t-D  Pi  p2  •••  PK 


aK,K  1 


(11. 3) 

(11. 4) 


Then, 


a0,K+l  "PK+1  a0,K 


(II. 5) 


r,K+l 


=  -p^  ,  a  +  a  ,  0  <  r  <  K 

*K+1  r,K  r-l,K 


(II. 6) 


aK+l,K+l  1 


(II. 7) 


Proof :  Equations  (II. 5)  and  (II. 7)  are  trivial  to  verify.  Equation  (II. 6)  is 
obtained  by  noting  that,  apart  from  the  sign,  R+1  is  a  sum  of  terms 

•  P 


p.  p. 

U  *12 


where  the  ii  ,  i2 ,  . . . ,  i„  are  all  possible  combinations 
XK+l-r  K+1_r 


of  K+l  -  r  subscripts  chosen  from  the  set  {l,  2,  ...,  K+ll .  This  is  equivalent 
to  summing  all  terms  p^  p^  ...  where  ix  ,  i2 ,  are  all 

possible  combinations  of  K+l  -  r  subscripts  from  the  set  |l,  2,  ...,  K|  and  then 

adding  this  to  the  sum  of  all  terms  p.  p. 

ll  12 

all  possible  combinations  of  K-r  subscripts  from  the  set  f  1,  2,  ...,  K|  multiplied 

by  p  .  Apart  from  the  sign  these  two  sums  are  none  other  than  a  and 
K+l  r-i,K 

a  multiplied  by  p  respectively.  Since  a  has  a  sign  term, 

X  y  K  K+ 1  X ,  IV 

K-r 

(-1)  ,  a  negative  sign  must  be  included  in  equation  (II. 6). 


p .  where  it  ,  i2  ,  . . . ,  i  are 
xK-r 
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Lemma  2:  Suppose  that  P.  j  =  0,  1,  _ _  K  are  defined  in  terms  of  the 

J  > 

px ,  p2 ,  . . . ,  p  as  above.  Let  pi  ,  p2 ,  . . . ,  p  be  zeroes  of  the  polynomial  of 
K  K. 

degree  K  in  X,  see  equation  (H . 1)  below 


K-l 


fK(X)  =  X  +  Vl,KX  +  +  ai,K  X  a0. 


Then  the  quantities  P.  can  be  expressed  in  terms  of  the  polynomial  coefficients 

J 


as 


Pj,K 


=  (-1) 


K-j 


( K-r  \  . 

ar,K  \j  -t)’  J 


0,  1,  2,  . . . ,  K  (II. 8) 


r=0 


Proof:  The  proof  is  by  induction.  For  K  =  2  the  P4  v  are  defined  by 

1  >  ^ 


Po,2  =  Pi  P2 

Pi  ,2  *  Pi  +  p2  -  2  Pi  p2 

P2  ,2  =  1  -  Pi  -  P2  +  Pi  P2 


and  the  polynomial  coefficients  a  are  defined  by 

r,K 


ao  j2 

=  Pi  P2 

al  ,2 

=  -(Pi  +  Pa) 

3-2  3  2 

=  1 

Simple  algebraic  rearrangement  yields 


?0  j2  =  &o  >2 

Pi  ,2  -  -(2ao,2  +  ai,2) 

P2  ,2  =  ao  ,2  +  ai  ,2  +  a2  ,2  • 


which  are  none 
To  complete 
K  =  H  and  show 


other  than  equations  (I I. 8)  for  K  =  2. 

the  inductive  proof  assume  that  equations  (II. 8)  are  true  for 
they  are  true  for  K  =  H+l.  Now 


Pj,H+l  ~  Pj,H  PH+1  +  Pj-1,H  (1  ‘  PH+1} 


PH+1  ^Pj,H  ~  +  Pj-1, 


H 


(II. 9) 


where  1  <  j  <  H.  Using  equations  (II. 8)  for  1  <  j  <  H, 
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P  -P 
J,H  j  -1,H 


j 


r=0 


a 


r,H 


(-1) 


H+l-j 


3-1 


r=0 


=  (-1) 


H-j 


r,H 


H+l-r 

j-T 


r=0 


since 


Substituting  this  in  equation  (II. 9) 


Pj  >H+1 


PH+1  *-  1') 


H-j 


r,H 


H+l-r 

j-r 


(-1) 


H+l-j 


j-1 


r,H 


H-r 
j-l-r. 


r=0 


r=0 


Using  the  result  of  Lemma  1 


Pj  jH+1 


(-1) 


H+l-j 


r  ,H+1 


H+l-r 

j-r 


r=0 


for  1  <  j  <  H.  For  j  =  0 


For  j  =  H+l, 


0,H+1 


=  Pi  P2  ...  p. 


■'  ^H+l 


,  i^H+1 

(-1)  a 


0,H+1 


PH+1,H+1  PH,H  ^  PH+1^ 


H 


ar,H  ^  "  Ph+1^ 


r=0 


H  H 


r=0  r=0 
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Using  Lemma  1  again 


'H+1,H+1 


H 


_PH+1  a0,H  +  2^  ar,H+l  +  aH,H 
r=l 


H+l 


ar,H+l  Since  aH,H  =  1  £°r  a11  H- 


r=0 


Lemma  3:  If  a  and  P.  are  defined  as  in  Lemma  2  then 
*  r ,  K.  j  ,  K. 


r,K 


(-1) 


K-r 


K-j 
j,K  Vr-j 


3=0 


(11.10) 


Proof :  Notice  that  equations  (II. 8)  can  be  manipulated  to  give  a  recursive  set 

of  equations  for  a  in  terms  of  P.  „  as 
r ,  K.  3  ,K 


0,K 


(-1)K  P 


0,K 


(II. 11) 


r-1 


r,K 


(_l)K_r  P 


r,K 


a. 

J,K  \r~3 


(11.12) 


j=0 


The  proof  of  equation  (II. 10)  is  inductive.  Obviously  equation  (II. 10)  is  true 
for  r  =  0.  It  is  shown  that  if  equation  (II.  10)  is  true  for  r  <v  then  it  is 

true  for  r  =  v.  Thus,  from  equation  (11.12)  for  v  >  1 


v- 1 


v,K 


(-1)K"V  P 


v,K 


,  K-j 

Ej  >K  l^v-j 


(11.13) 


j=0 


In  the  summation  of  equation  (11.13)  the  following  substitution  can  be  made 


j,K 


(-DK-j  Vp  fK-m' 

Lj 
m=0 


which  results  in 
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v- 1 


v,K 


C-1JK-V  P 


v,K 


(-D 


k-j  v  d  ( K-mY K-j 


m. 


,K  Vj-m/Vv-j 


j=0 


m=0 


=  (-D 


K-v 


v- 1 


v,K 


(-D 


j=0 


v-j  \  p  ( K-m  V  K-j 

m,K  \  j-m/ \v-j 

m=0 


Using  the  fact  that 


v  3 


j  =0  i=0 


’ij  =  l_j  bJ  i__j  8Ji 
j=0  i=j 


v ,  K 


=  c-i) 


K-v 


v,K 


V-  1  V- 1 


Pj  >K 


j=0  m=j 


(-D 


v-m  [  K-j 
m-j 


K-m 

v-m 


Now 


v- 1 


(-1) 


v-mf  K-j 
m-j 


K-m 

v-m 


v- 1 


K-j 

v-j 


m=J 


m=J 


(-1) 


v-mf  v-j 
m-j 


since 


K-j  K-m 
m-j  /\v-m 


(K-j)!  (K-m)! 


(K-m) !  (m-j)  !  (K-v) !  (v-m) ! 


(K-j)! 


(v-j)  ! 


(v-j)!  (K-v)!  (v-m)!  (m-j)! 


K-j  Vv-j 
v-j  /V  m-j 


If  it  is  shown  that 


v—  1 


(-1) 


v-m  (  v-j  \ 


m-j 


=  -1 


m=J 


the  proof  is  complete.  In  this  summation  put  m  =  b  +  j  and  the  summation  becomes 
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v-j-1 


b=0 


Using  the  binomial  theorem,  (a  +  b)n 


a^n  ^  >  with  a 
j=0 


-1,  b  =  1  gives 


n-1 


j=0 


2(v-j)-l  _ 


Thus,  the  summation  is  equal  to  (-1) 


-1. 
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TABLE  1.  STATE  PROBABILITIES  FOR  FIXED  PRI  EMITTERS 


Gate  time  (r)  limits 

W  =  0,1,2. . . 

Possible  number  of 
pulses 

States  of  Model 

Model  state 
probabilities 

Wz, 

W 

0 

„  _  CW+l)z  -  r 
P°  -  2 

r  -  Wz 

(W+l)  z 

W+l 

1 

P>  =  z 

TABLE  2.  STATE  PROBABILITIES  FOR  TWO  POSITION  STAGGER  EMITTERS,  z2  >  zx 


Case 

Gate  time  (r)  limits 

W  =  0,1,2... 

Possible 

number 

of 

pulses 

States 

of 

model 

Model 

state  probabilities 

W(zx+z2), 

2W 

0 

Po  = 

(2W+1) (zj  +z2 )  -  2 t 

Zl+Z2 

A 

W  ( Zj  +  z2  )  +  z, 

2W+1 

1 

Pi  = 

2r  -  2W(Zi+z2) 

Zj+Z2 

W(Zj  +z2  )  +  Zj  , 

2  W 

0 

Po  = 

W(zi  +z2 )  +  z2  -  r 

Zl  +z2 

B 

2W+1 

1 

Pi  = 

2zj 

Zl  +z2 

W(zi+z2)  +  z2 

2W+2 

9 

p2  = 

r  -  zj  -  W(zi  +z2 ) 

Zi  +z2 

C 

W(zx+z2)  +  z2  , 

2W+1 

0 

Po  = 

(2W+2) (zi +Z2 )  -  2 t 

Zl+Z2- 

(W+l)Cz1+z2) 

2W+2 

l 

Pi  = 

2r  -  (2W+1)  (Zi  +z2 ) 

Zl  +z2 
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TABLE  3(a).  STATE  PROBABILITIES  FOR  THREE  POSITION  STAGGER  EMITTERS, 

Z3  >  Z2  >  Zl  ,  Zl  +  Z2  >  Z3 


Gate  time  (r)  limits 

W  =  0,1,2. . . 

Possible 

number 

of 

pulses 

W(z,  +z2+z3) , 

3W 

W(Zi  +z2+z3)  +  zi 

3W+ 1 

W(zi  +z2  +z3)  +  zi  , 

3W 

3W+1 

W(  zi+z2+z3)  +  z2 

3W+2 

W(zi+z2+z3)+z2  , 

3W 

3W+1 

W(Zi  +z2  +z3  )  +  z3 

3W+2 

W(z,  +z2+z3)  +  z3  , 

3W+1 

W(Zi  +Z2  +Z3  )  +Zi  +z2 

3W+2 

W(Zi+Z2+Z3)+Zi+Z2  , 

3W+ 1 

3W+2 

W(Zi  +Z2+Z3)+Zi  +z3 

3W+3 

W(Zi+Z2+Z3)+Zi  +z3. 

3W+1 

3W+2 

W(Z!+Z2+Z3)+Z2+Z3 

3W+3 

W(zi  +z2  +z3  )+z2  +z3  , 

3W+2 

(W+l)  (Zi+Z2+Z3) 

3W+3 

States 

of 

model 


Model  state  probabilities 


(3W+1)  (zj  +  z2  +z3  ) -3r 
Zl  +z2+z3 

3t-3W(z!+z2+z3) 

Zj  +z2  +z3 

z3 +z2 -2r +2W(zi +z2 +z3  ) 

Z l  +Z2+Z3 

r  +2zi  -W ( Zi  +z2  +z3  ) 

”  Zl  +z2+z3 

T  -  Zl  -W(Zi  +z2+z3) 

Zl  +  z2+z3 


Z3 -T+W(Zi  +z2+z3) 

Zl  +Z2  +z3 

2  (Zi  +z2  )  -r  +W(zj  +z2  +z3  ) 

Zl  +  Z2  +  Z  3 

2r -zi -z2 -2W(zj +z2 +z3 ) 
Zi  +z2  +z3 


(3W+2)  (zi  +z2  +  z3  ) -3t 
Zl  +z2+z3 

3r  -  (3W+1)  ( z i  +z2  +z3  ) 

Zi  +z2+z3 

2  (z3  -t  )  +  Zi  +z2  +2W(zi  +z2  +z3  ) 
Zl  +z2  +z3 

r-z3  +  zi  +z2 -W(zi  +z2+z3  ) 

Zi+Z2+Z3 

T-Zi  -Z2-W(Z!+Z2+Z3) 

Zl  +Z2  +Z3 

Z2+Z3-T+W(Zi+Z2+Z3) 


Z]  +z2+z3 

3zi  +z2+z3  -r+W(zi  +z2  +  z3) 


Zl  +z2  +  z3 

2(r-zj  )-z2  -z3  -2W(zj  +z2+z3) 


Zl  +z2  +z3 


(3W+3)  (zi +z2 +z3  ) -3r 

Zi+Z2+Z3 

3r  -  (3W+2)  (zi +z2 +z3  ) 
Zl  +z2  +z3 
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TABLE  3(b) .  STATE  PROBABILITIES  FOR  THREE  POSITION  STAGGER  EMITTERS, 

Z3  >  Z2  >  Zi  ,  Zi  +  Z2  <  Z3 


Case 

Gate  time  (r)  limits 
W  =  0,1,2. . . 

Possible 

number 

of 

pulses 

States 

of 

model 

Model  state  probabilities 

1 

W(z1+z2+z3), 

3W 

0 

Po 

(3W+1)  (zi  +z2  +z3  ) -3r 

Zj+Z2+Z3 

■ 

W(zi+z2+z3)+zj 

3W+1 

1 

Px 

3t-3W(z1+z2+z3) 

Zi+Z2+Z3 

W(zi+z2+z3)+zi, 

3W 

0 

Po 

Z3  +  z2-2r+2W(zi+z2+z3) 

Zi  +z2+z3 

B 

3W+1 

1 

Pi 

_  r+2zi -W(z3 +z2+z3 ) 

Zi+Z2+Z3 

W(z1+z2+z3)+z2 

3W+2 

2 

P2 

_  T-Z!  -W(zx  +z2  +z3) 

Zj+Z2+Z3 

W(Zj  +z2  +z3)+z2  , 

3W 

0 

Po 

z3-r+W(z!+z2+z3) 

Zi  +z2  +z3 

C 

3W+1 

1 

Pi 

2  (zj  +z2  )  -r  +W  (zi  +z2  +z3  ) 

Zj+Z2+Z3 

W(Z!+Z2+Z3)+Zi+Z2 

3W+2 

2 

P2 

=  2t-zj  -z2  -2W(zi  +z2+z3  ) 

Zl  +Z2  +Z3 

W(Zi+Z2+Z3)+Z!+Z2  , 

3W 

0 

Po 

Z3-T+W(Zj  +z2+z3) 

Zl  +z2+z3 

D 

3W+1 

1 

Pi 

Z!+Z2 

Zx  +z2  +z3 

3W+2 

2 

P2 

Zl  +z2 

Zl  +z2  +z3 

W  (zj  +z2  +z3  ) +z3 

3W+3 

3 

P3 

_  T-Zi -Z2-W(Zi+Z2+Z3) 

Zi+Z2+Z3 

W(z3  +z2+z3)+z3  , 

3W+1 

0 

Po 

2  (Z3  -T  )  +zx  +z2  +2W  (Zl  +z2  +z3  ) 

Zi+Z2+Z3 

E 

3W+2 

1 

Pi 

=  r-z3+zi+z2-W(zi+z2+z3) 

Zi+Z2+Z3 

W(z3  +z2+z3)+z3  +z3 

- — - 

3W+3 

2 

P2 

_  T-Zi-Z2-W(Zi+Z2+Z3) 

Zl  +Z2+Z3 
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TABLE  3(b) (CONTD.) . 


Case 

Gate  time  (r)  limits 

W  =  0,1,2. . . 

Possible 

number 

of 

pulses 

States 

of 

model 

Model  state  probabilities 

W(zj  +Z2  +  Z3 )"*"Z1  +z3  , 

3W+1 

0 

Po 

Z2+Z3 -T+W(Z! +Z2+Z3) 

Zj  +  z2  +  z3 

F 

3W+2 

1 

Pi 

3zi  +Z2+Z3  -r+W(zx  +Z2+Z3) 

Zj  +Z2  +Z3 

W(Zi  +  Z2+Z3)+Z2+Z3 

3W+3 

2 

P2 

2(r-zi  )-z2 -z3-2W(zi+Z2+z3) 

Z1+Z2+Z3 

W(Zi+Z2+Z3)+Z2+Z3  , 

3W+2 

0 

Po 

(3W+3)  (z!  +z2  +Z3 )  -3t 

Z1+Z2+Z3 

G 

(W+l)  (z1+z2  +  z3) 

3W+3 

1 

Pi 

3r  -  (3W+2)  (zj  +z2  +z3  ) 

Zj  +z2  +z3 
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TABLE  5.  STATE  PROBABILITIES  FOR  EMITTERS  WITH  UNIFORM  PRI 
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TABLE  6.  PROBABILITIES  FOR  PRI  TYPE  RECOGNITION  AND  PARAMETER  ESTIMATION 


Type  of 

PRI 

variation 

r  <  a 

T  =  r J  =  (a+b)/ 2 

t  =  2Ti  =  a+b 

Fixed 

Po  =  (Z-T)/z 

Po  =  0 

Po  =  0 

PRI 

p2  =  0 

P2  =  0 

p2  =  0 

Zi  +  z2  -2t 

Z2  -Zj 

Two 

position 

Po  =  7  +7 

Zj  +z2 

P°  2  (zi  +z2  ) 

Po  =  0 

stagger 

P2  =  0 

P2  =  Po 

p2  =  0 

a+b-2r 

b-a 

b-a 

Po  "  a+b 

Po  4 (a+b) 

Po  3 (a+b) 

Uniform 

p2  =  0 

P2  =  Po 

P2  =  Po 

a+b-2r 

b-a 

n  4  (b-a) 

Po  ~  a+b 

Po  ~  7r  (a+b) 

Po  jt2  (a+b) 

Sinusoidal 

P2  =  0 

P2  =  Po 

P2  =  Po 
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Figure  1 
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*  denotes  the  first  pulse  after  start  of  gate  interval  t 


Figure  1.  Graphical  derivation  of  state  probability  formulae 
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Figure  3 


Filter  1 


Pulse 

detector 


Filter  2 


Pulse 

detector 


Figure  3.  Detecting  overlapping  pulses 


Pulse 

counter 
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Figure  4 


Emitter  PRI's  are  z1  =  103  Ms  and  z2  =  223  Ms  (nominally) 


Figure  4.  Effect  of  r  and  (Tr )  on  occurrence  of  complex  zeroes 
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Figure  5 


Emitter  PRI's  are  z1  =  103  Ms  and  z2  =  223  Ms  (nominally) 


Figure  5.  Effect  of  r  and  (Tr)  on  PRI  estimate  error 
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Figure  6 


Emitter  PRI's  are  z1  =  103  Ms  and  z2  =  223  Ms  (nominally) 


Figure  6.  Disciminant  of  f2 (X)  for  simulated  environment 
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PRI  range  is  200  -  5000  Ms 

Number  of  runs  is  1000  -  Gate  time,  t  =  100  Ms 
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Figure  7.  Performance  of  nonadaptive  PRI  estimates 
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Figure  8.  Performance  of  adaptive  estimate  of  equations  (49),  (50) 
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Figure  9.  Performance  of  adaptive  estimate  of  equations  (49),  (51) 
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Figure  10.  Probability  of  error  for  adaptive  classification  of  emitters 
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Figure  11.  Performance  of  adaptive  parameter  estimates 
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Figure  12.  Model  verification  results 
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Figure  13.  Model  verification  for  greater  PRI  range 
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Figure  14.  Model  verification  for  two  position  stagger  PRI 
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